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ABSTRACT 



A sample of 20 powerful extended radio galaxies with redshifts between zero 
and two were used to determine constraints on global cosmological parameters. 
Data for six radio sources were obtained from the VLA archive, analyzed, and 
combined with the sample of 14 radio galaxies used previously by Guerra & Daly 
to determine cosmological parameters. The new results are consistent with our 
previous results, and indicate that the current value of the mean mass density 
of the universe is significantly less than the critical value. A universe with Q m 
in matter of unity is ruled out at 99.0% confidence, and the best fitting values 
of Vt m in matter are O.lOl^io an d — 0.25^25 assuming zero space curvature and 
zero cosmological constant, respectively. Note that identical results obtain when 
the low redshift bin, which includes Cygnus A, is excluded; these results are 
independent of whether the radio source Cygnus A is included. The method 
does not rely on a zero-redshift normalization. 

The radio properties of each source are also used to determine the density 
of the gas in the vicinity of the source, and the beam power of the source. The 
six new radio sources have physical characteristics similar to those found for the 
original 14 sources. The density of the gas around these radio sources is typical 
of gas in present day clusters of galaxies. The beam powers are typically about 
10 45 erg s _1 . 

Subject headings: cosmology: observations — galaxies: active, evolution, jets - 
radio continuum: galaxies 
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1. Introduction 

The future and ultimate fate of the universe can be predicted given a knowledge of 
the recent expansion history of the universe (assuming the universe is homogeneous and 
isotropic on scales greater than the current horizon size). This recent expansion history 
can be probed by studying the coordinate distance to sources at redshifts of one or two; 
the coordinate distance is equivalent to the luminosity distance or angular size distance 
multiplied by factors of (1 + z). The advantage of determining cosmological parameters 
using the coordinate distance (luminosity distance or angular size distance) is that this 
distance depends on global, or average, cosmological parameters. It is independent of the 
way the matter is distributed spatially, of the power spectrum of density fluctuations, of 
whether the matter is biased relative to the light, and of the form or nature of the dark 
matter (assuming that the universe is homogeneous and isotropic on large scales). 

It was shown in 1994 that powerful double-lobed radio galaxies provide a modified 
standard yardstick that can be used to determine global cosmological parameters (Daly 
1994, 1995), much like supernovae can be used as modified standard candles. The method 
was applied and discussed in detail by Guerra & Daly (1996, 1998), Guerra (1997), and 
Daly, Guerra, & Wan (1998) who found that the data strongly favor a low density universe; 
a universe with Vt m = 1 was ruled out at 97.5 % confidence. 

It was shown by Daly (1994, 1995) that the radio properties of these sources could be 
used not only to study global cosmological parameters, but also to determine the ambient 
gas density, beam power, Mach number of lobe advance, and ambient gas temperature of the 
sources and their environments. The characteristics of the sources and their environments 
are presented and discussed in a series of papers (Wellman & Daly 1996a,b; Wan, Daly, & 
Wellman 1996; Daly 1996; Wellman, Daly, & Wan 1997a,b; Wan & Daly 1998a,b; Wan, 
Daly, & Guerra 2000). 
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Radio maps of six powerful double-lobed radio galaxies were extracted from the Very 
Large Array (VLA) archives at the National Radio Astronomy Observatory (NRAO), and 
analyzed in detail. New results on global cosmological parameters, ambient gas densities, 
and beam powers are presented here. 

The expanded sample is described in §2. The new results are presented in §3. The 
implications of the results are discussed in §4. 

2. Expanded Sample 

Each powerful extended radio galaxy (also known as a "classical double") has a 
characteristic size, D*, that predicts the lobe-lobe separation at the end of the source 
lifetime (Daly 1994, 1995; Guerra & Daly 1998). The parameters needed to compute the 
characteristic size are the lobe propagation velocity, vl, the lobe width, ol, and the lobe 
magnetic field strength, B L . These three parameters can be determined using radio maps 
with arc-second resolution at multiple frequencies, such as those produced with the VLA or 
MERLIN (Multi-Element Radio Linked Interferometer Network). Multiple- frequency data 
are needed to use the theory of spectral aging to estimate the lobe propagation velocity 
(e.g., Myers & Spangler 1985). In addition, these maps must have the necessary angular 
resolution and dynamic range to image sufficient portions of the radio bridges. 

Two published data sets, Leahy, Muxlow, & Stephens (1989) and Liu, Pooley, & Riley 
(1992), have radio maps of powerful extended radio galaxies at multiple frequencies which 
are sufficient to compute all three parameters used in determining D*. These data were 
used to compute for 14 radio galaxies (Guerra & Daly 1996, 1998; Guerra 1997; Daly, 
Guerra, & Wan 1998). Current efforts to expand the data set are underway, and include 
searches through the VLA archive. The VLA archive search has already yielded the desired 
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data for six radio galaxies, and new results including these six sources are presented here. 

Data were selected from observations of powerful extended radio galaxies from the 3CR 
sample (Bennett 1962) on the basis of the observation frequency and array configuration 
used. An observation in the VLA archive was a candidate if the lobe-lobe angular size of 
the source was 10 to 40 times the implied beam size, and four hours separated the first 
and last scans. These observations should resolve the source sufficiently and have enough 
■uf-coverage to image the bridges. From candidate observations at both L and C band, we 
have successfully imaged data from six sources. The VLA archive data sets used here are 
listed in Table 1. 

Radio imaging was performed using the NRAO AIVS software package. The uv data 
needed minimal editing, and initial calibration was performed in the standard manner using 
AIVS. As diagnostics, initial images were made using the AIVS task IMAGR both without 
and with the CLEAN deconvolution algorithm. The final radio maps were produced with the 
SCMAP task, which performs self-calibration along with the IMAGR and CLEAN tasks. 

Parameters used in AIVS for imaging, deconvolution, and calibration are chosen to 
produce radio maps for a given source that can easily be used for spectral aging analysis. 
This is particularly important where observations at different frequencies are produced by 
different observers. Although some of these radio maps exist, different choices of parameters 
(such as the restoring beam) are often used to produce these radio maps. Thus, reducing 
the raw uv data insures that the data analysis can be performed consistently between radio 
maps. 
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2.1. Radio Maps from Archive Data and Their Use 

Radio intensity maps produced from VLA archive data are presented here (Figures 
1-6). The 4.872 GHz map of 3C 244.1 was created with a restoring beam identical to the 
synthesized beam of the 1.411 GHz map of 3C 244.1, and the 4.885 GHz map of 3C 194 
was created with a restoring beam identical to the synthesized beam of the 1.465 GHz map 
of 3C 194. Thus, radio maps at both frequencies for a particular radio galaxy have similar 
angular resolutions. Radio maps created with three of these archive data sets have appeared 
in the literature. A 1.411 GHz map of 3C 244.1 appeared in Leahy & Williams (1994), and 
4.88 GHz maps of 3C 244.1 and 3C 325 appeared in Fernini, Burns, & Perley (1997). 

The computation of a^, B L , and vl were performed here in the same manner as 
Wellman, Daly, & Wan (1997a,b). The deconvolved lobe width, a L) is measured lO/i -1 
kpc from the hot spot toward the host galaxy. The lobe magnetic field strength, B L) is 
computed from the deconvolved surface brightness and bridge width measured 10/i" 1 kpc 
from the hot spot toward the host galaxy. The lobe propagation velocity, Vl, is computed 
on the basis of spectral aging along the imaged portions of the bridge, and the magnetic 
field used in spectral aging is computed using values measured lOh' 1 kpc and 25h~ 1 kpc 
from the hotspot. The errors on each of these quantities are discussed in §5 of WDW97b, 
§5.2 of WDW97a, and the appendix of this paper. 

All parameters presented here are computed with b = 0.25, which means magnetic 
fields are computed to be 0.25 times the minimum energy values, and without an a — z 
correction which refers to a correction related to the observed correlation between spectral 
index and redshift (see Wellman 1997; Wellman, Daly, Wan 1997a,b; Guerra 1998; Guerra 
& Daly 1998 for details). It was found by Guerra & Daly (1998) that constraints on 
cosmological parameters did not depend on these choices, and very similar results are 
obtained independent of the value of b and of whether ana-z correction is applied. 
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2.2. Theory 

The small dispersion in the average size of 3CR radio galaxies at a given redshift 
suggests that, at a given redshift, all of the sources have a very similar average size. This 
size may be estimated by the average size of the full population of powerful extended radio 
galaxies at that redshift, < D >, or by the average size of a given source at that redshift, 
-D*. If the total time a source produces powerful jets, with beam power Lj, that powers 
the growth and radio emission of the source is then the average source size will be 
-D* = v L t*, assuming the rate of growth of the source, vl, is roughly constant over the 
source lifetime. The source velocity vl, estimated via synchrotron and inverse Compton 
aging techniques, increases with redshift (Leahy, Muxlow, & Stephens 1989; Liu, Pooley, & 
Riley 1992; Daly 1994, 1995; WDW97b), and there is no indication that the rate of growth 
of the sources decreases with redshift. However, the average size of the full population 
decreases monotonically with redshift for redshifts greater than about 0.5 (see Table 3 in 
this paper, and Figure 1 in Guerra & Daly 1998). 

This means that, for sources of this type, must decrease with redshift. For the 
sources studied here, it was shown by WDW97a that the radio power of a given source is 
roughly constant over the lifetime of that source, making it very unlikely that the radio 
power of a given source decreases with time, causing sources at higher redshift to fall below 
the radio flux limit of the survey when they are smaller. 

Thus, since < D > decreases with redshift for z > 0.5, and D* = vl £* with vl 
either increasing with redshift or independent of redshift, the data require that t* decrease 
with redshift. The rate of growth of the source, vj_,, depends on one parameter that is 
intrinsic to the AGN, the beam power Lj, and on two parameters that are extrinsic to 
the AGN, the ambient gas density n a and the cross sectional area of the radio lobe oc a\: 
vl oc {Lj n~ l a^ 2 )^ 1 /^. The total time for which the AGN produces a powerful outflow, t*, 
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must depend on properties intrinsic to the black hole, and the only intrinsic parameter that 
appears in vl is Lj. Thus, we write oc LJ 13 ^ 3 , which defines the one model parameter (3. 
The other intrinsic property of the black hole that might enter is the total energy available 
to power the outflow, E 1 *, but since E* — Lj t*, any dependence on E* is absorbed into the 
relation between t* and Lj. 

A power law relation between the total time the AGN produces powerful jets and the 
beam power of the jets is not unexpected. It means that there is a power law relation 
between the beam power Lj and the total energy available initially it is assumed that 
the beam power is roughly constant over the lifetime of the source, but the value of Lj is 
set by the initial energy available E*. This is reminiscent of the power law relation for 
main sequence stars between luminosity and lifetime, or between total energy available and 
luminosity for main sequence stars, and of other power law relations that arise so frequently 
in astrophysics. For example, for a value of (3 of 2, the beam power of the jet is related to 
the total energy available by the relation Lj oc El, which is similar to the relation between 
luminosity and mass for main sequence stars: L oc M 3 ' 5 . Note, that for main sequence stars, 
the luminosity is determined by the total energy available initially, and remains roughly 
constant over the lifetime of the main sequence star. For a value of (3 of 1.5, the relation 
between beam power and total energy available is Lj oc E%. 

The excellent fits obtained by comparing < D > with assuming that t* oc Lj' 3 ^ 3 
supports this choice for the parameterization of t*. As shown in Figure 7, if the zero-redshift 
bin containing Cygnus A is excluded from the analysis, and the best fit value of (3 and the 
constant C rat i (defined below) are determined, we can predict the value of for Cygnus 
A, and it matches the value of for Cygnus A to very high accuracy. This is quite 
extraordinary if we consider the large drop in < D > that occurs from the first redshift bin 
to the zero-redshift bin, which is mirrored by the large drop in of Cygnus A (see Figures 
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8 and 9). 

The theory requires that track < D > with redshift, thus it requires 

< D > / D* = C ra tio , 

where C ratio is a constant, independent of redshift. We do NOT require that < D > /D* — 1, 
but we require that this ratio is a constant, C rat i . Thus, the use of this equation for 
cosmology is independent of any factor that might be used to normalize the quantity D*. 

The ratio < D > / D* depends on cosmological parameters and the model parameter 
(3. The way that this ratio depends on coordinate distance is given by equation (6) of 
Guerra & Daly (1998), and the full dependence of D* on observable parameters is given in 
Appendix A of this paper. 

The confidence with which cosmological parameters and the model parameter can 
be constrained depends on the uncertainty of the ratio < D > / D*. This is obtained 
by combining the uncertainty of < D > and the uncertainty of -D* in quadrature. The 
uncertainty of < D > is given by the dispersion in radio source size at a given redshift. The 
errors in < D > range from 10 % to 35 % (see Table 3). This dominates the uncertainty 
of the ratio < D > /D*. The uncertainty of -D* depends on the uncertainty of each of the 
quantities used to determine D*, and is discussed in detail in Appendix A of this paper. 



2.3. Application of the Theory 

In terms of quantities that can be estimated from radio observations, Lj oc vl B\ a 2 L 
so oc vl LJ 13 ^ 3 oc t>^ 1 ~ /3 ^ 3 ' ) (B L a L )~ 2 ^^ 3 for perfectly symmetric lobes and bridges. The 
three parameters, a^, B^, and vl are computed for the bridge on each side of the six 
radio galaxies obtained from the VLA archive, with the exception of one bridge in 3C 324 
which was not imaged along its length sufficiently The characteristic core-lobe size, r*, was 
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computed for each bridge using using the equation described above and in Guerra & Daly 



The characteristic (lobe-lobe) size, D*, is taken to be the sum of both r* (or in the case of 
3C 324, = 2r*), and is normalized so that D* of Cygnus A (3C 405), a very low redshift 
source in our sample, is equal to the average size of the full population of powerful classical 
double radio galaxies at very low redshift; this does not impact the determination of the 
model parameter, j3, or constraints on cosmological parameters, in any way (see §2.2). 
Constraints on cosmological parameters determined using this method are independent of 
the normalization of D*, as described in §2.2. As discussed in §3.1, the best fit value of the 
one model parameter (5 is determined simultaneously with the best fit values for the two 
cosmological parameters that enter, Q m and Qa, the normalized values of the current values 
of the mean mass density and the cosmological constant (see Guerra & Daly 1998; Daly, 
Guerra, & Wan 1998). 

Table 2 presents the six new D t values in the last column, assuming the best fit value of 
P = 1.75 (see §3.1). Source name and redshift are listed in the first two columns. The third 
column lists the redshift bin corresponding to the assignments in Guerra & Daly (1998) 
and Table 3 below. The lobe-lobe angular size of the source is listed in the fourth column, 
the fifth and sixth columns list the core- lobe characteristic sizes, r*, and the characteristic 
source size is listed in column seven. 



(1998), and Daly (1994, 1995): 




(1) 
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3. Updated Results 

3.1. Constraints on Cosmological Parameters 

The subsample of sources with estimates of the characteristic size, D*, has been 
increased from 14 to 20, as discussed above in §2. Only sources with physical sizes, defined 
as the projected separation between the radio hot spots, greater than 20 hr x kpc can 
be used to determine a characteristic size. This is because smaller sources are typically 
not sufficiently resolved so that the radio data are useful, and the radio lobes of smaller 
sources are interacting with the interstellar medium of the host galaxy rather than the 
intergalactic/intracluster medium. It was decided that this same criterion should be applied 
to the larger comparison sample of powerful 3CR radio galaxies. Thus, the sample of radio 
galaxies used to determine the redshift evolution of the physical size has been reduced from 
82 to 70; twelve radio galaxies were cut because the physical separation between their lobes 
was less than 20 h~ x kpc. This has a rather small impact on the actual means and standard 
deviations of the parent population, as can be seen by comparing Table 3 of this paper 
with Table 1 of Guerra & Daly (1998). The average lobe-lobe separations as a function 
of redshift are listed in Table 3 for three example choices of cosmological parameters 
(matter-dominated, curvature-dominated, and spatially flat with non-zero cosmological 
constant). 

To solve simultaneously for the model parameter (5 and the cosmological parameters 
VL m and ^a, the ratio of D* for each source to (D), the average lobe-lobe size of the 
parent population in the corresponding redshift bin, is fit to a constant, independent of 
redshift: < D > /D* — C rat i 01 where the constant C rat i is allowed to float when the 
best fit parameters are determined. The value of the constant is a free parameter, so the 
normalization of D* does not affect the fits in any way. 
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Figure 7 illustrates the cosmological dependence of (D)/D* on the coordinate distance 
(a Q r). As described below, the best fit value of (3 obtained here is (3 = 1.75; for this 
value of (3, (D)/D* is proportional to (a r) L6 . Thus ((D) / D :¥ )(a r)~ ls is independent of 
cosmological parameters. The data can be compared to several different sets of cosmological 
parameters on a single figure by plotting ((D) / ' D*)(a r)~ l - % for each data point and 
comparing this with (a r) -1 - 6 curves obtained for different sets of cosmological parameters, 
as is shown in Figure 7. Note that, for the fits shown, the lowest redshift bin, which includes 
Cygnus A, was excluded from the analysis. Even so, the lines all go directly through this 
point (denoted by a star) indicating the predictive power of this model. 

The hypothesis is that, for the correct choice of cosmological parameters, (D)/D* = 
Cratio = constant, so that the values of ((D) / D*)(a r)~ ls for all 20 radio galaxies should 
follow a curve that, at each z, tracks the curve (a r)~ ls obtained for that particular choice 
of cosmological parameters. Figure 7, shows ((D) / D 1r )(a r)~ 1 - & for the the six new points 
and 14 original points as a function of z (the six new points are denoted by open squares). 
Also drawn on this figure are the best-fit curves of (a r)~ L6 obtained for specified values of 
cosmological parameters and excluding Cygnus A from the fits. In this figure, all of the 
curves pass through Cygnus A, though this is not required when we actually 
solve for best fitting cosmological parameters. Including Cygnus A has a negligible 
effect on the normalization of these curves and a small effect on cosmological constraints. It 
is clear that curves obtained for a low density universe describe the data points quite well, 
and the curve describing a universe with Q m = 1 does not follow the data points. 

For all 20 sources, the chi-squared for fitting the ratio to a constant is computed for 
relevant values of (3, f2 m , and f^. It is found that the best fit value of f3 is f3 — 1.75 ± 0.25. 
This result is insensitive to the choice of Q m and Q^, and there appears to be no significant 
covariance between (3 and cosmological parameters (see Figures 10a and 10b). 
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The confidence contours in the Vt m - Vt^ plane are shown in Figures 11 and 12. The 
probability associated with a given range of Q m and fl\ independent of j3 is shown in Figure 
11 (referred to as two-dimensional confidence intervals). In Figure 12, the projection of a 
confidence interval onto either axis (Q m or Q\) indicates the probability associated with 
a given range of that one parameter, independent of all other parameter choices (referred 
to as one-dimensional confidence intervals). Both figures illustrate how this method and 
the data are most consistent with a low density universe; Vt m < 0.15 with 68% confidence, 

~ 0.5 with 90% confidence, and Vl m < 1.0 with 99% confidence. The constraints on f2 A 
are not as strong, and values of Q\ from zero to unity are consistent with the data. 

The best fit value of (3 = 1.75 ± 0.25 is consistent with the previous estimates of 
P ~ 1.5 ± 0.5 (Daly 1994), and (3 ~ 2.1 ± 0.6 (Guerra & Daly 1996, 1998), but with 
significantly reduced uncertainties. Similarly, the constraints on cosmological parameters 
are consistent with previous estimates, (Daly 1994; Guerra & Daly 1996, 1998; Guerra 
1997; Daly, Guerra, & Wan 1998) but with smaller error bars. It is apparent in Figures 
11 and 12 that these data and method strongly favor a low density universe; a universe 
where Vt m = 1 is ruled out with 99.0 % confidence independent of ^a an d (5 . This will 
be discussed in more detail in §4. To further illustrate the insensitivity of cosmological 
constraints to including Cygnus A, Figures 13 and 14 show the confidence intervals when 
Cygnus A is excluded from the fits. 

Figures 13 and 14 also illustrate that that data are most consistent with a low density 
universe. Fits excluding Cygnus A rule out an Vt m = 1 universe at 97.5% confidence. 
These results are quite similar to those obtained when Cygnus A is included. This further 
illustrates that the method does not rely upon a low-redshift normalization. 
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3.2. Ambient Gas Densities and Beam Powers 

Daly (1994, 1995, 2000), following the work of Perley & Taylor (1991), Carilli et. al 
(1991), and Rawlings & Saunders (1991), showed that the radio properties of a powerful 
extended radio source could be used to estimate the beam power, Lj, and density of the 
gas in the vicinity of the source, n a . The method is described in detail by Wellman (1997), 
Wellman, Daly, & Wan (1997a), Wan (1998), Wan, Daly, & Guerra (1998), and Daly (2000). 
Values for the ambient gas density and beam power for the 14 sources in the original sample 
are described in these papers. The basic equations are: 

Lj oc a| P L v L , (2) 

and 

n a oc ; (3) 

where Pl is the lobe pressure (see equation A2), and the normalizations are given in the 
references listed above. The values obtained for the six new radio galaxies in our sample 
are listed in Table 4. Also listed in Table 4 are all the input parameters used to compute 
Lj, n a , and D*. 



4. Discussion 

A parent population of 70 powerful extended classical double radio galaxies with 
redshifts between zero and two was used to define the evolution of the mean or characteristic 
size of these sources as a function of redshift; these sources are all Type 1 FRII sources 
(Leahy & Williams 1984), also referred to as FRIIb sources (Daly 2000). An independent 
estimate of the mean or characteristic size of a given source was possible for a subset of 
20 of these radio galaxies for which extensive multiple frequency radio data was available. 
Requiring that the two measures of the mean source size have the same redshift behavior 
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allows a simultaneous determination of three parameters: the one model parameter (3, 
and the two cosmological parameters Q m and Q A (assuming that the only significant 
cosmo logical parameters today are the mean mass density, a cosmological constant, and 
space curvature). 

The method was applied to this data set, and interesting new constraints are presented. 
It is found that the model parameter is very tightly constrained to be (3 = 1.75 ± 0.25 (see 
Figure 10), consistent with previous estimates, and that this model parameter is independent 
of cosmological parameters. For a value of (3 = 1.75, the characteristic source size is 

oc (B L a L )~ 7 / 6 v L 5 / 12 , which indicates that it is necessary to have multiple-frequency 
radio data in order to estimate -D*, owing to its vl dependence. 

The data strongly favor a low density universe; a universe with Q m = 1 is ruled out 
with 99% confidence, independent of the value of or (3. Either space curvature or a 
cosmological constant, or both, are allowed. The main conclusion is that Q m is low, but, 
at this point, the method and data do not allow a discrimination between whether space 
curvature or a cosmological constant is important at the present epoch. 

It is interesting to note that the lowest reduced chi-squared obtained is 0.96 for 
VL m ~ —0.25 and il\ ~ 0. This value is slightly greater and closer to unity than the 
minimum reduced chi-squared obtained by Guerra & Daly (1998), which indicates that this 
sample of 20 sources has a reasonable distribution around any model predictions. This 
convergence to unity with increasing sample size suggests that this method and its statistics 
are reliable. 

The best fit for cosmological parameters in the physically relevant half-plane of Q m > 
is VL m = and = 0.45 with a reduced chi-squared of 0.98 (also close to unity). However, 
Figures 11 and 12 clearly show that our results are still consistent with VL\ = 0. 



-16- 



The radio data also allow a determination of the density of the ambient gas in the 
vicinity of each radio source, and the beam power of each source; the values of these 
quantities are presented. The sources lie in high-density gaseous environments like those 
found in low-redshift clusters of galaxies. Typical beam powers for the sources are 
~ 10 45 erg s _1 . 

Special thanks go to Rick Perley and Miller Goss for their aid in extracting data from 
the VLA archive. The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated Universities, Inc. 
The authors would also like to thank Katherine Blundell, Chris Carilli, Paddy Leahy, Wil 
van Breugel, and Dave Wilkinson for helpful discussions. Special thanks are extended to the 
referee for providing valuable feedback and suggestions that led to significant improvement 
of the paper. This work was supported in part by the U.S. National Science Foundation 
and the College of Liberal Arts and Sciences at Rowan University. 



A. Error Analysis 

Errors on and errors on < D > both contribute to the confidence contours 
relevant for the cosmological parameters determined here. The errors on < D > are purely 
statistical, and arise from the dispersion in source size of the parent population, and at the 
present time dominate the total uncertainty of the quantity < D > /D*. The errors on D* 
depend on the uncertainties of quantities used to determine D*; the way that the error on 
D* is computed is presented here. 

The quantity 

r, = v L U<x v L LjW oc v L (P L v L a\) ~ m oc v\~^ a?'* P?' 3 . (Al) 
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The lobe pressure 

Pl oc (± b- 15 + b^ Bl iw (A2) 

where B min is the minimum energy magnetic field, B min oc SV 7 aT^ 1 , and is determined 
using parameters obtained 10 h~ x kpc behind the hotspot: the radio surface brightness at 
this location is S^io and the lobe half- width at this location is cll,w- The radio surface 
brightness and lobe half-width 25 h^ 1 kpc behind the hotspot are denoted 5^25 and ai,25 
respectively. The lobe propagation velocity vl is estimated using a standard synchrotron 
and inverse Compton aging model in which the time, t, required for the source to grow a 
size Ax, is used to estimate the lobe propagating velocity vl = Ax/t. The time 

B 1/2 

t oc - Bav (A3) 

4 /2 (Bl + B MB ) 

(see, for example, Wan & Daly 1998, §3) where v T is the break frequency, B av is the average 



magnetic field in the radio bridge, taken to be B av = V-B10-B25 (WDW97b), B w is the 
magnetic field strength 10 hr x kpc behind the hotspot, -B25 is the magnetic field strength 
25 hr l kpc behind the hotspot, and B MB is the term that describes inverse Compton 
cooling of relativistic electrons by scattering with microwave background photons, and is 
obtained by equating the energy density of the microwave background at a given redshift 
with B MB /(8n). 

In terms of these parameters, 

/ A \ -13/3 

r. oc (1 + J B W W a L ^ 3 Ax^^ vf~ m (A4) 

v / R 2 , R 2 \ 1-/3/3 R -|(l-/3/3) R -i(l-/?/3) 
x \ D av -°MSj -°10 -°25 

Given that the magnetic field strength is parameterized by B = b B m i n at any given 
location, we obtain the following expression for r*, which is relevant for a determination of 
the error: 

n oc (V L5 + 6 2 j b~^ 1+ V Ax^-W VT ^-m bJW 3 -V a L)10 -^ (/3 - 1/7) (A5) 
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o -09/6+1/14) 
A *-V,10 



(b 2 S v , 10 



2/7 C 2/7 _ -2/7 _ -2/7 , d 2 
^,25 a L,10 °L,25 + -£>MB 



(1-/3/3) 



The total error on (57%, /r* is found by taking the partial derivative of r\ with respect 
to a given variable, for each variable, adding these terms in quadrature, and taking the 
square- root. We obtain 



f ^1 
V r* 



5r, 
r* 



Ax 



(3/3) 
(1-/3/3) 



5r 



B 2 s 



-4(1-/3/3) + (1-/3/3) 



Ax y 

-B10-B25 



a L,10 



i(l/7-/3)-|(l-/3/3) . 2 



-B10-B25 + ^Mfl' 
-B10-B25 



25 



S//,l() 

2/3 



-W + 1/14) + |d - /J/3) ( . ^" B 'l , 

I V-D10-D25 + H MB 



B25 

£ql,io 

a L,10 
Sv,10 



1_ & 3.5 ^ X 

' -^l + Z?) + (1-/3/3) 



2^x0-025 



(A6) 
(A7) 
(A8) 
(A9) 
(A10) 
(AH) 



• VI 3 • l> : '>- r 'J 2" \B 10 B 25 + B MB 

The final uncertainty of divided by r*, 5r*/r* is obtained by adding each of the 
terms listed above in quadrature, and taking the square root. As mentioned above, this 
is typically much less than the uncertainty of < D >, so the confidence level with which 
cosmological parameters are determined is primarily controlled by the uncertainty in < D > 
(see Table 3). 

The uncertainties on each of the quantities listed above are described in section 5 
of WDW97b, and section 5.2 of WDW97a. The uncertainty in v T is typically (10 to 20) 
%. The uncertainty on Ax is typically (6 to 20) %. The uncertainty on B 2 s is typically 
6 %. The uncertainty on a Ljl0 is typically (6 to 20)%; this term (et. A9) is generally the 
largest contribution to the error in r*. The uncertainty on SV,io is typically 20 %, and the 
uncertainty on b is about 15 % for b = 0.25, and is about 10 % for b = 1 (see Wellman, 
Daly,& Wan 1997b, section 7.4). 
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For a given value of /3, equations A6 through All may be written as a number times 
the fractional error on the relevant quantity. Some of the equations involve the fraction 
/ = (^ b w bII+b mb ^ ) ' wn i cn i s always less than or equal to one. To get a sense of which 
terms contribute to the total error in r*, equations A6 through All are rewritten assuming 
a value of (3 of 1.75. 



<h ' A 0.10 + 0.42/) (A14) 



^) = (-0.80- 0.12/) (^^1 (A15) 

(^) =(-0.36 + 0.12/) f^^l (A16) 

(S H -<- u » + "fl(T) < A17 » 

= (-0.514- 0.83/) ff) (A18) 



V r * / b=0.25 V ^ 

The largest contribution to the uncertainty in r* is usually due to the uncertainty in 
the bridge width Ql,io- 
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Table 1. VLA Archive Data Sets 









Listed 


Frequency 


Array 




source 


r>ana 


Program ID 


Observer 




Config. 


iJate 


3C 244.1 


L 


POOL 


Pooley, G. 


1.411 


B 


09/17/82 




C 


AF213 


Fernini, I. 


4.872 


B 


12/23/91 


3C 337 


L 


AR123 


Rudnick, L. 


1.452 


A 


02/21/85 




C 


AR123 


Rudnick, L. 


4.885 


B 


06/01/85 


3C 325 


L 


AV153 


van Breugel, W. 


1.465 


A 


12/05/88 




C 


AF213 


Fernini, I. 


4.885 


B 


12/23/91 


3C 194 


L 


AV164 


van Breugel, W. 


1.465 


A 


05/11/90 




C 


AV164 


van Breugel, W. 


4.885 


A 


05/11/90 


3C 324 


L 


AR123 


Rudnick, L. 


1.452 


A 


02/21/85 




C 


AR123 


Rudnick, L. 


4.885 


B 


06/01/85 


3C 437 


L 


AV164 


van Breugel, W. 


1.465 


A 


05/11/90 




C 


AV164 


van Breugel, W. 


4.885 


B 


04/21/89 
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Table 2. Radio Galaxies with D* Presented Here 









e 




a 




Source 


z 


Bin 


(arcsec) 


(h- 1 kpc) 


(h- 1 kpc) 


(h- 1 kpc) 


3C 244.1 


0.43 


2 


50.8 


151 ± 16 


135 ± 16 


286 ± 22 


3C 337 


0.63 


3 


43.5 


150 ± 24 


65 ±7 


214 ±25 


3C 325 


0.86 


3 


15.8 


164 ± 52 


66 ± 14 


230 ± 54 


3C 194 


1.19 


4 


14.2 


105 ± 22 


92 ± 15 


197 ±27 


3C 324 b 


1.21 


5 


10.2 


75 ± 18 




149 ± 35 


3C 437 


1.48 


5 


36.7 


55 ±6 


48 ±5 


103 ±8 



a Computed assuming (3 = 1.75, Ho = 0.1, = 0, b = 0.25 and not including the a — z 

correction. 

for only one bridge. 



Table 3. The Average Lobe-Lobe Sizes for Powerful 3CR Radio Galaxies. 



(D) (h' 1 kpc) 



Bin z Range Sources fi G = 1.0, fi A = 0.0 fi G = 0.1, fi A = 0.0 fi G = 0.1, fi A = 0.9 



1 


0.0-0.3 


3 


66±14 


68±14 


72±13 


2 


0.3-0.6 


13 


202±45 


224±50 


259±57 


3 


0.6-0.9 


23 


148±17 


174±20 


209±24 


4 


0.9-1.2 


16 


107±24 


133±29 


165±36 


5 


1.2-1.6 


9 


91±32 


122±43 


152±53 


6 


1.6-2.0 


6 


66±19 


92±26 


114±33 
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Table 4. Summary of New Source Properties 



oource 






b 

a L 


D C 
-DlO 


D d 

-D25 


v L e 


n a 


log Lf 


244.1 


0.43 


101 


5.5 ± 0.3 


3.5 ± 0.2 


2.5 ± 0.2 


1.1 ± 0.2 


1.0 ± 0.4 


44.35 ± 0.10 






91.6 


4.7 ± 0.4 


5.3 ± 0.3 


3.2 ± 0.2 


1.7 ± 0.3 


1.0 ± 0.4 


44.75 ± 0.11 


337 


0.63 


125 


6.2 ± 0.5 


3.5 ± 0.2 


2.3 ± 0.1 


1.5 ± 0.5 


0.6 ± 0.4 


44.57 ± 0.15 






69.8 


10.5 ± 0.4 


3.8 ± 0.2 


2.8 ± 0.2 


1.1 ± 0.2 


1.2 ± 0.4 


44.99 ± 0.10 


325 


0.86 


49.0 


3.3 ± 1.1 


6.8 ± 0.9 


4.0 ± 0.3 


2.0 ± 0.5 


1.2 ± 0.6 


44.73 ± 0.23 






34.7 


4.4 ±0.8 


12.7 ± 1.0 


4.6 ±0.3 


3.0 ±0.8 


1.9 ± 1.0 


45.70 ±0.17 


194 


1.19 


31.6 


4.5 ± 1.0 


7.1 ±0.6 


5.0 ±0.3 


1.8 ±0.3 


1.6 ±0.6 


45.00 ±0.16 






50.0 


5.6 ±0.8 


6.7 ±0.5 


5.6 ±0.4 


2.0 ±0.3 


1.2 ±0.4 


45.16 ±0.13 


324 


1.21 


22.1 


4.5 ±0.9 


10.3 ±0.9 




2.1 ±0.7 


2.4 ± 1.5 


45.38 ±0.19 


437 


1.48 


106 


11.4 ±0.8 


6.9 ±0.4 


4.2 ±0.2 


4.6 ±0.8 


0.2 ±0.1 


46.18 ±0.10 






101 


10.5 ±0.8 


9.4 ±0.6 


6.7 ±0.4 


6.0 ± 1.0 


0.2 ±0.1 


46.49 ±0.10 



a Core-hotspot distance in hr x kpc. 

b Lobe radius, 10 h" 1 kpc behind hot spot, in hr x kpc. 

c Minimum energy magnetic field, 10 hr x kpc behind hot spot, in /i 2//7 10~ 5 G. 
d Minimum energy magnetic field, 25 hr l kpc behind hot spot, in /i 2//7 10 _5 G. 
c Lobe advance speed, in 10 _2 c. 

f Ambient gas density, in 10~ 3 /i n cm -3 , where n = 12/7 for B » Bmb and n = 20/7 for 
B « B M b- 

g Logarithm of the luminosity in directed kinetic energy, in h~ 2 erg s -1 . 



-26- 



(a) 

PLot file version 2 created 21-MAY-2000 16:47:37 
CONT:3C244.1 IPOL 1411.250 MHZ POOL244 1.ICLN.1 



58 30 45 
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30 
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00 






Jq 




29 45 






n 




30 




o 


\ J 

I I I 


i i 



10 30 24 22 20 18 16 14 

RIGHT ASCENSION 
Cont peak flux = 1.4806E+00 JY/BEAM 
Levs = 1 .000E-02 * (-1 , 1 , 2, 4, 8, 1 6, 32, 64, 
128) 




Fig. 1. — Radio maps of 3C 244.1 from archive data in (a) L band and (b) C band. 
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PLot file version 2 created 21-MAY-2000 16:15:50 
CONT:3C337 IPOL 1477.400 MHZ 337LA.ICLN.1 
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RIGHT ASCENSION (B1950) 
Cont peak flux = 9.1149E-01 JY/BEAM 
Levs = 2.500E-03 * (-1, 1, 2, 4, 8, 1 6, 32, 64, 
128, 256) 
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PLot file version 2 created 21-MAY-2000 16:19:21 
CONT:3C337 IPOL 4860.100 MHZ 337CB.ICLN.1 



44 25 55 




16 27 21 



19 18 17 

RIGHT ASCENSION (B1950) 
Cont peak flux = 3.6257E-01 JY/BEAM 
Levs = 5.000E-04 * (-1 , 1 , 2, 4, 8, 1 6, 32, 64, 
128, 256, 512) 



Fig. 2. — Radio maps of 3C 337 from archive data in (a) L band and (b) C band. 
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(a) 

PLot file version 2 created 19-MAY-2000 14:46:13 
CONT:3C325 IPOL 1489.900 MHZ 325LA.ICLN.1 

I I I I I I 




(b) 

PLot file version 2 created 19-MAY-2000 14:39:17 
CONT:3C325 IPOL 4860.100 MHZ 325CB.ICLN.1 

I I I I I I 




Fig. 3. — Radio maps of 3C 325 from archive data in (a) L band and (b) C band. 
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(a) 

PLot file version 2 created 19-MAY-2000 14:04:28 
CONT: 0806+426 IPOL 1489.900 MHZ 194LA.ICLN.1 




O 



08 06 39.0 38.5 38.0 37.5 37.0 

RIGHT ASCENSION (B1950) 
Cont peak flux = 6.5124E-01 JY/BEAM 
Levs = 2.000E-03 * (-1 , 1 , 2, 4, 8, 1 6, 32, 64, 
128, 256) 



(b) 

PLot file version 5 created 19-MAY-2000 14:21 :09 
CONT: 0806+426 IPOL 4860.100 MHZ 194CA.ICLN.1 

n i i r 




Fig. 4. — Radio maps of 3C 194 from archive data in (a) L band and (b) C band. 
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PLot file version 2 created 19-MAY-2000 14:29:41 
CONT: 3C324 IPOL 1477.400 MHZ 324LA.ICLN.1 
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(b) 



PLot file version 2 created 19-MAY-2000 14:18:02 
CONT:3C324 IPOL 4860.100 MHZ 324CB.ICLN.1 
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Fig. 5. — Radio maps of 3C 324 from archive data in (a) L band and (b) C band. 
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PLot file version 2 created 21-MAY-2000 16:32:22 
CONT: 2145+151 IPOL 1489.900 MHZ 437LA.ICLN.1 
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(b) 



PLot file version 2 created 21-MAY-2000 16:24:53 
CONT:3C437 IPOL 4860.100 MHZ 437CB.ICLN.1 
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Fig. 6. — Radio maps of 3C 437 from archive data in (a) L band and (b) C band. 
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Cygnus A 
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z 

Fig. 7. — The quantity ((D) / 'D*)(a r) _1 - 6 , computed assuming (3 = 1.75 and b = 0.25. This 
measured quantity is independent Q m and Q\. Open squares denote the six new data points 
from the VLA archive and Cygnus A is denoted by the star symbol. For different choices 
of Q m and Q\, the best-fit redshift evolution of (a r) -1 ' 6 excluding Cygnus A is shown. 
Note that the fits are not normalized to Cygnus A, but are consistent with Cygnus A. 
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Fig. 8. — The average lobe-lobe size (D) of powerful extended 3CR radio galaxies computed 
assuming Q m = 0.1 and fi^ = 0. 
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Fig. 9. — The quantity computed assuming (5 = 1.75, b = 0.25, Vt m = 0.1, and = 0. 
Open squares denote the six new data points from the VLA archive and Cygnus A is denoted 
by the star symbol. 
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Fig. 10. — The 68% and 90% confidence intervals for (3 and Q m , where (a) ^a=0 (no 
cosmological constant)and (b) il\ = 1 — Q m (spatially flat universe). 
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Fig. 11. — The 68% and 90% confidence intervals for ranges of both fl m and independent 
of (3. (Two-dimensional) 
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Fig. 12. — The projections of the 68%, 95%, and 99% confidence intervals onto either 
axis (Q m or indicates the probability associated with the range in that one parameter, 
independent of all other parameter choices. (One-dimensional) 



